Mult if r act al statistics of Lagrangian velocity and acceleration in turbulence 
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The statistical properties of velocity and acceleration fields along the trajectories of fluid particles 
transported by a fully developed turbulent flow are investigated by means of high resolution direct 
numerical simulations. We present results for Lagrangian velocity structure functions, the acceler- 
ation probability density function and the acceleration variance conditioned on the instantaneous 
velocity. These are compared with predictions of the multifractal formalism and its merits and 
limitations are discussed. 
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Understanding the Lagrangian statistics of particles 
advected by a turbulent velocity field, u(x,t), is impor- 
tant both for its theoretical implications [1] and for ap- 
plications, such as the development of phenomenological 
and stochastic models for turbulent mixing [2] . Recently, 
several authors have attempted to describe Lagrangian 
statistics such as acceleration by constructing models 
based on equilibrium statistics (see e.g. [3, 4, 5], critically 
reviewed in [6]). In this letter we show how the multi- 
fractal formalism offers an alternative approach which is 
rooted in the phenomenology of turbulence. Here, we de- 
rive the Lagrangian statistics from the Eulerian statistics 
without introducing ad hoc hypotheses. 

In order to obtain an accurate description of the par- 
ticle statistics it is necessary to measure the positions, 
X(t), and velocities, v(t) = X(i) = u(X(i),t), of the 
particles with very high resolution, ranging from frac- 
tions of the Kolmogorov timescale, t v , to multiples of the 
Lagrangian integral time scale, Tl- The ratio of these 
timescales, Tt,/t 71 , gives an estimate of the micro-scale 
Reynolds number, R\, which may easily reach values 
of order 10 3 in laboratory experiments. Despite recent 
advances in experimental techniques for measuring La- 
grangian turbulent statistics [7, 8, 9], direct numerical 
simulations (DNS) still offer higher accuracy albeit at a 
slightly lower Reynolds number [10, 11, 12, 13]. In this 
letter we are concerned with single particle statistics, that 
is, the statistics of velocity and acceleration fluctuations 
along individual particle trajectories. Here, we analyse 
Lagrangian data obtained from a recent DNS of homo- 
geneous isotropic turbulence [14] which was performed 
on 512 3 and 1024 3 cubic lattices with Reynolds numbers 
up to R\ ~ 280. The Navier-Stokes equations were in- 
tegrated using fully de-aliased pseudo-spectral methods 
for a total time T « Tl- Millions of Lagrangian particles 
(passive tracers) were released into the flow once a statis- 
tically stationary velocity field had been obtained. The 



positions and velocities of the particles were stored at a 
sampling rate of 0.07t^. The Lagrangian velocity was 
calculated using linear interpolation. Acceleration was 
calculated both by following the particle and by direct 
computation from all three forces acting on the parti- 
cle - the pressure gradients, viscous forces and the large 
scale forcing. The two measurements were found to be in 
very good agreement. The flow was forced by keeping the 
total energy constant in the first two wavenumber shells 
(for more details see [14]). 

It is well known that Lagrangian velocity increments, 
S T v = v(t+T)—v(t), are quasi-Gaussian for time lags r of 
order Tr, but become increasingly intermittent at higher 
frequencies [8]. The resulting acceleration statistics ex- 
hibit some of the most extreme fluctuations of any known 
quantity, with accelerations, a(t), up to 80 times its root 
mean square value, a rmsi possible [7]. The most natu- 
ral way to quantify such phenomena is via probability 
density functions (pdfs) of the Lagrangian velocity in- 
crement, V(S T v), and acceleration, V(a). The frequency 
of extreme events is reflected in the size of the tails of 
the pdfs and thus in the high order moments. These can 
be analysed with the aid of Lagrangian velocity struc- 
ture functions S p (t) — {(S T v) p ) where S T v characterises 
the magnitude of a component of the velocity increment. 
Since the flow here is isotropic the choice of component 
is immaterial. 

It has long been recognised that Eulerian velocity fluc- 
tuations in the inertial subrange exhibit anomalous scal- 
ing: ({S r u)P) = ({u(x + r) - u{x)Y) ~ [15], where 
r is the spatial separation. On the basis of simple phe- 
nomenological arguments, we may expect the Lagrangian 
velocity fluctuations to exhibit a power law behaviour for 
time scales within the inertial subrange too. We may 
therefore assume that S p (t) ~ t^ l ^ with t v «C t <C Tl. 
Anomalous scaling is often interpreted as the result of the 
intermittent nature of the energy cascade. Among the 
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simplest stochastic models able to reproduce both quali- 
tatively and quantitatively such intcrmittency are those 
based on the multifractal formalism. This has been suc- 
cessfully used to explain Eulcrian statistics such as struc- 
ture functions [15, 16, 17] and velocity gradients [18, 19] 
and Lagrangian statistics such as the acceleration covari- 
ance [20] and the velocity statistics [21, 22]. The aim 
of this letter is to compare predictions of the multifrac- 
tal formalism for Lagrangian velocity structure functions, 
the acceleration pdf and the acceleration variance condi- 
tioned on the instantaneous velocity with those obtained 
from the DNS data. The Lagrangian multifractal pre- 
dictions are derived from the multifractal formalism in 
the Eulerian reference frame without any additional free 
parameters. 

In the multifractal formalism the global scale invari- 
ance of Kolmogorov's theory (K41) becomes a local scale 
invariance. Namely, the turbulent flow is assumed to 
possess a range of scaling exponents I = (h m i n , h max ). 
For each h S I there is a set Sh € M 3 of fractal dimen- 
sion D{h) such that, in the limit of small r: S r u(x) ~ 
Uo(r/Lo) h ^ for x G Sh- Here mq is the large scale 
fluctuating velocity and L is the integral length scale. 
For small values of uq we are in the laminar part of the 
flow for which a multifractal description is not appro- 
priate. From this local scaling law, the scaling proper- 
ties of the Eulerian structure function can easily be de- 
rived by integrating over all possible h [15]: ((5 r u) p ) ~ 
(ug) j I dh(r/L ) hp+3 - D{h '>. The factor (r/Lo) 3- ^ is 
the probability of being within a distance of order r of 
the set Sh of dimension D(h). A saddle point approxima- 
tion in the limit rdo then gives the scaling exponents 



( E (p) = -mf(hp + 3-D(h)). 
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(1) 



If the Eulcrian scaling exponents are known, D(h) can 
be calculated from the inverse of the Legendre transfor- 
mation (1). Among possible empirical formulae for the 
scaling exponents, (e(p), we choose the one of She and 
Leveque [23] . Using this it can be shown that 



D(h) = l+p*{h) 
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(2) 



where p*(h) = (3/ ln(2/3)) ln[(l - 9h)/(6 hi(2/3))] is the 
value of p which minimises the inverse of (1). 

The velocity fluctuations along a particle trajectory 
may be considered as the superposition of different con- 
tributions from eddies of all sizes. In a time lag r the 
contributions from eddies smaller than a given scale, r, 
are uncorrelated and one may then write 5 T v ~ 5 r u. We 
assume that r and r are linked by the typical eddy turn 
over time at the given spatial scale, r r ~ r/S r u. There- 
fore, in the multifractal terminology, 



S T v ~ 5 r u 



Lo i- 
r ~ — r 



(3) 



The presence of fluctuating eddy turn over times is the 
only additional complication introduced by the multifrac- 
tal formalism in the Lagrangian reference frame. Using 
(3) we can now derive a prediction for the Lagrangian 
velocity structure function [22]: 

hp+3-D(h) 

T 

ihei 

where the factor (t/Tl)^- ^^ 1 -^ is the probability 
of observing an exponent h in a time lag r. The expo- 
nents (l(p) then follow from a saddle point approxima- 
tion in the limit r -C Tj,: 

' hp + 3 - D(hY 



S P (r) ~ <«g) / dh 
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(4) 
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FIG. 1: ESS plot of Lagrangian velocity structure function 
S p (t) versus S^r). Symbols refer to the DNS data for 
p — 8, 6, 4 from top to bottom. Lines have slopes Ci(p)/C^(2) 
given by the multifractal prediction (4). In the inset we show 
the local slopes of the DNS data and the multifractal predic- 
tions versus the time lag, r. 

In Fig. (1) the results for S p (t) calculated from the 
DNS are presented. Although the scaling in a log-log plot 
is reasonable, a more detailed inspection of logarithmic 
local slopes displays a deterioration of scaling properties 
at small times. This is due to the presence of a strong 
bottleneck for time lags, r s [t v , 10t ti ]. This has been ex- 
plained in terms of trapping events inside vortical struc- 
tures [14], a dynamical effect which may strongly affect 
scaling properties and which a simple multifractal model 
cannot capture. For this reason, scaling properties are 
recovered only using Extended Self Similarity (ESS) [24] 
and for large time lags, r > 10T n . In this interval a sat- 
isfactory agreement with the multifractal prediction (4) 
is observed, namely Cl(4)/Cl(2) = 1.71; Cl(6)/Cl(2) = 
2.16;Cl(8)/Cl(2) = 2.72. 

Similar phcnomenological arguments can be used to 
derive predictions for the acceleration statistics. The ac- 
celeration at the smallest scales is defined by 

a = -—. (5) 
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As the Kolmogorov scale, 77, fluctuates in the multifrac- 
tal formalism [15]: rj{h,v a ) ~ (vL l [ ! l /vo) 1 ^ 1+h ^ , so does 
the Kolmogorov time scale, T v (h,Vo)- Using (3) and (5) 
evaluated at 77, we get, for a given h and u , 



a(h,v ) ~ v i+*v^ +h L 



_ _ 3h 

1 + h j 1 + h 



(6) 



The pdf of the acceleration can be derived by integrat- 
ing (6) over all h and v , weighted with their respective 
probabilities, (r„(/i, «o)/7lM) (3 ~- D(/i))/(1 ~' i) &ndT(v ). 
The large scale velocity pdf is reasonably approximated 
by a Gaussian [15]: V(vq) = cxp(— Vq/2(7 2 )/ ^/2tw^, , 
where cr 2 = (vq). Integration over w gives: 

spy/ \ f AU fe ' 5 + D(fe) 7 - 2fe - 2P " 1 > T D(h) + h-3 -1 

V(a) ~ / d/ia 3 ^ 3 CT « x 



exp -- 



2(l + fe) 2(l-2fe) „ 

a 3 v 3 Ln 



(7) 



From (7) we can derive the Reynolds number de- 
pendence of the acceleration moments [20, 25]. For 
example, in the limit of large R\ the second or- 
der moment is given by (a 2 ) cx R\ where x — 
swp h {2(D(h) -Ah- l)/(l + /i)). Thus, we find that 
X = 1-14 which differs slightly from the K41 scaling, 
X K41 = 1 (see [25, 26] for a discussion on departures from 
K41 scalings in the context of acceleration statistics). In 
order to compare the DNS data with the multifractal 
prediction we normalise the acceleration by the rms ac- 
celeration, a a = (a 2 ) 1 / 2 . In terms of the dimensionless 
acceleration, a~aja a , (7) becomes 



V{a) 



(h-5 + D(h)) 



R v x {h) exp 
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1_2(1 + M z , h s 

-a 3 r 



dh, 
(8) 

where y(h) = x{h-5+D(h))/6+2(2D(h) + 2h-7)/3&nd 
z(h) = x(l + h)/S + 4(2ft - l)/3. We note that (8) may 
show an unphysical divergence for a w for many mul- 
tifractal models of D(h). For example, with D{h) given 
by (2) we cannot normalise P(a) for h < h c w 0.16. 
This shortcoming is unimportant for two reasons. First, 
as already stated, the multifractal formalism cannot be 
trusted for small velocity and acceleration increments be- 
cause it is based on arguments valid only to within a con- 
stant of order one. Thus, it is not suited for predicting 
precise functional forms for the core of the pdf. Sec- 
ond, values of h < h c correspond to very intense velocity 
fluctuations which have never been accurately tested in 
experiments or by DNS. The precise functional form of 
D{h) for those values of h is therefore unknown. Thus, 
we restrict h to be in the range h c < h < h max . For 
hmnx we take the value of h which satisfies D'(h) = 0, 
that is h max w 0.38. Values of h > h max affect only the 
peak of the velocity distribution which we have already 
excluded from our discussion. We also restrict |a| to lie 
in the range [a m i n , 00) with a m ; n = O(l). 
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FIG. 2: Log-linear plot of the acceleration pdf. The crosses 
are the DNS data, the solid line is the multifractal predic- 
tion and the dashed line is the K41 prediction. The DNS 
statistics were calculated along the trajectories of two million 
particles amounting to 3.6 x 10 9 events in total. The statis- 
tical uncertainty in the pdf was quantified by assuming that 
fluctuations grow like the square root of the number of events. 
Inset: a 4 V(a) for the DNS data (crosses) and the multifractal 
prediction. 



In Fig. (2) we compare the acceleration pdf com- 
puted from the DNS data with the multifractal pre- 
diction (8). The large number of Lagrangian particles 
used in the DNS (see [14] for details) allows us to de- 
tect events up to 80 a a . The accuracy of the statis- 
tics is improved by averaging over the total duration of 
the simulation and all directions since the flow is sta- 
tionary and isotropic at small scales. Also shown in 
Fig. (2) is the K41 prediction for the acceleration pdf 
V Kil {&) - a- 5 / 9 i?^ 1/2 exp(-a 8 / 9 /2) which can be re- 
covered from (8) with h = 1/3, D(h) = 3 and X Kil = 1- 
As is evident from Fig. (2), the multifractal prediction 
(8) captures the shape of the acceleration pdf much bet- 
ter than the K41 prediction. What is remarkable is that 
(8) agrees with the DNS data well into the tails of the 
distribution - from the order of one standard deviation 
a a up to order 70cr a . This result is obtained with D(h) 
given by (2). We emphasise that the only degree of free- 
dom in our formulation of V(a) is the minimum value 
of the acceleration, a m i n , here taken to be 1.5. In the 
inset of Fig. (2) we make a more stringent test of the 
multifractal prediction (8) by plotting a 4 V(a) and which 
is seen to agree well with the DNS data. 

From (6) it is also possible to derive a prediction for 
the acceleration moments conditioned on the local - in- 
stantaneous - velocity field vo: (a n \vo). For example, for 
the conditional acceleration variance we get 
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l+4h-D(h) 3 + D(h) £>(h)-6h-3 

dhv r+s „ i+" l r+s _ 



In the limit v <5£ 1, a saddle point approximation gives 
(a 2 \v ) cx l/av (3+D(h))/(i+h) where a = . nf ^ ((1 + Ah _ 
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D(h))/(1 + h)) and h is the value of h which minimises 
the exponent of v. Thus, we find that (a 2 |«o) oc Vq' 57 . 



ramctcr in our formulation is the value of a m i n which docs 
not have a sensitive effect on the results. 



1000 




The simulations were performed on the IBM-SP4 of 
Cineca (Bologna, Italy). We thank in particular G. Er- 
bacci and C. Cavazzoni for resource allocation and pre- 
cious technical assistance. We acknowledge support from 
the EU under the contracts HPRN-CT-2002-00300 and 
HPRN-CT-2000-0162. 



FIG. 3: Log-log plot of the conditional acceleration variance. 
The crosses are the DNS data, the solid line is the multifrac- 
tal prediction and the dashed line is the prediction of [25]. 
The DNS data is the absolute acceleration conditioned on the 
absolute velocity. Statistical uncertainty was estimated by 
dividing the samples into five sub-ensembles. Inset: the con- 
ditional acceleration variance scaled by the multifractal pre- 
diction, (a 2 \v) /v 4 ' 57 , and the prediction of [25], (a 2 \v) /v 6 , 
both normalised by a 2 .. 

In Fig. 3 we plot (a 2 \v) 1 normalised by the accelera- 
tion variance, versus v 2 /°fl- The relatively large error 
that can be seen in the DNS conditional acceleration 
statistics for large values of v 2 /a 2 reflects the rarity of 
these events. However, in agreement with [25] a clear 
trend is evident that, for large velocities, the accelera- 
tion magnitude depends strongly on the magnitude of 
the velocity. (The vector acceleration and velocity can 
easily be shown to be uncorrelated for stationary tur- 
bulence [25].) Also shown in Fig. 3 are the multifractal 
prediction and the prediction of [25] based on a dimen- 
sional argument pertaining to the vorticity, namely that 
(ct 2 |t>) oc v e . (This can be recovered from the multifrac- 
tal formalism with h — 0, corresponding to a singularity 
in S T v, and D(h) = 3.) Although statistical noise pre- 
vents us from making a convincing claim, the multifractal 
prediction appears to agree better with the DNS data. 

In conclusion, we have shown that the multifractal for- 
malism predicts a pdf for the unconditional acceleration 
which is in excellent agreement with the DNS data. In 
deriving the form of the pdf we have only assumed that 
the Lagrangian velocity increment is related to the Eu- 
lerian velocity increment by (3) and that the large scale 
fluctuating velocity is Gaussian. The only adjustable pa- 
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